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Under this NASA contract, the University of Alabama in Hunts- 
ville (UAH) was to provide the following: 

1. identification of critical flow problems related to 
propulsion systems and identify computational fluid 
dynamics (CFD) resources for application to these 
problems ; 

2. promotion of CFD technology om research centers to 
technology development centers and direct the applica- 
tion of CFD as a design tool; 

3. encourage industry and university participation in 
CFD/ propulsion research activities through their own 
internal and external funds; 

4. provide peer review to CFD/ propulsion programs; 

5. direct the verification and validation of CFD methodolo- 
gies as they are applied to propulsion problems. 

I terns 1~4 of the proposed UAH tasks above have been accomplished 
through the establishment of four Principle I n vest i ga tor Working 
Groups o Technology Teams under the Consortium for Computational 
Fluid Dynamics Application in Propulsion Technology (see Figures 
1 and 2) established by the CFD Branch at NASA/MSFC. Coordina- 
tion of these activities has been through the Consortium for 
Computational Fluid Dynamics (CCFD), established at and by UAH 
for this purpose. The four es t ab 1 i s hed technology teams are: 

- Turbine Stage Technology Team 

- Pump Stage Technology Team 

- Combustion-Driven Flow Technology Team 

~ Complex Flow Paths Technology Team 

The technical management of each of these teams is controlled by 
a staff member from the CFD Branch at NASA/MSFC and the members 
of each technical team are from the member organizations in the 
CCFD (see Figure 3). 

The published objectives and tasks of the CCFD are listed 
be 1 ow . 

Objectives : 


1. Focus CFD applications in propulsion 
a. ETC 

i , Direct baseline program towards improved a 
racy, stability and efficiency 
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b. CST l 

i. Stimulate CFD validation towards propulsion 
f 1 ows 

ii. Direct applications codes toward design tools 
and advanced hardware technology concepts 

2 - Identify national CFD propulsion requirements 

3. Stimulate a forum for government, industry, and univer- 
sity interactions 

4. Encourage industry to participate in CFD development 
with TRAD fund s 

5. Provide synergism in the CFD community 

6. Provide peer review of CFD programs 
T asks : 

1. Develop a plan to apply CFD to current and future pro- 
pulsion systems 

a. Identify and rank critic£tl flow problems related to 
propulsion systems 

b. Identify national CFD related resources 

c. Define high performance computing requirements to 
accomplish CFD for propulsion applications 

2. Direct CFD technology development to propulsion applied- 
t i o ns 

3. Assess and validate CFD applications In propulsion 
systems 

a. Develop evaluation criteria 

b. Define and implement benchmark validation 

c. Define and implement validation tests 

4. Direct the application of CFD design tools towards 
advanced hardware tec kino logy concepts 

5. Accelerate the transfer of CFD technology from universi- 
ties and research centers to industry and hardware 
development centers 

Three of the Technology Teams are well established and have been 
meeting for about one year; the fourth, the Complex Flow Paths 
Technology Team, is in start-up phase. These teams have been 
very useful in identifying current design issues ana promoting 
CFD solutions to these problems. The quarterly meetings have 
been useful in driving the progress of the work and coordinating 
the scheduling of results by various organizations. 

Items 1 and 5 of the proposed UAH tasks mentioned previously 
have been carried out by the principal investigator, one research 
associate and four graduate students at UAH. Three specific 



were accomplished to 


support these tasks. 


One effort i 
generation using 
carried out in 
explanation of how 
tional results are 


nvolved geometric modeling and 
a software tool called GENIE, 
collaboration with an engineer 
to use this software and a ser 
given in Appendix 1. 


numerical grid 
This work was 
at MSFC. The 
es of cornputa- 


Another effort to s 
sis and optimisation of 
this investigation and 
Appendix 2 . 


upport the CFD Branch has been the analy- 
a large CFD code, ROTOR. The results of 
the analysis of ROTOR are included in 


A third effort has been a 
chemistry model. The results of 
Appendix 3 . 


validation of an 
this analysis are 


equ i 1 i br ium 
included in 
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Recommendations 


It was agreed between UAH and MSFC that during the initial 
implementation of the four technology teams that the development 
of the Consortium Council (see Figures I and 2) be delayed. 

It is recommended that the implementation of the Council be 
considered for 1991. The Technology' Teams (Principal Investiga- 
tor Working Groups) will be well established by then and the 
overall propulsion program will have benefited from their ef- 
forts. This would also be the appropriate time to establish a 
communication network among government, industry and university 
consortium members. This network would include newsletters and 
brochu res . 
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APPENDIX 1 
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INTRODUCTION 


This report explains the activities that were carried out in 
order to generate a grid on the IRIS Silicon Graphics worksta- 
tion. The IRIS is a graphics workstation that allows the user to 
see a graphical re pres enta t i on of a solution for a particular 
problem. The grid generation tool that was used is called GENIE. 
GENIE is a 2-D and 3-D geometry modeling and grid generation tool 
that is available on the IRIS workstation . 
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IMPLEMENTATION 


^ To access the IRIS workstation, the user must input the 
login number and the password. The syntax for this interaction 


Pr omp t Input 

login: Enter login number and return 

password: Enter password and return 

The user gets a prompt on the screen that contains the name of 
the workstation with a number. The syntax of this prompt is: 

[ (deckard) or (rachel) ]. log name 8% ] 

In order to d i f f erent i ate the different workstations that are 
available, different names were given to them. The names of the 
ones that are available are deckard and rachel. For more infor- 
mation on the naming convention, the system administrator should 
be consu 1 t ed . 


To run GENIE or, the workstation, the user command is 
The IRIS workstation runs on the UNIX operating system 
means that all commands will be written in lower case 
because of UNIX sensitivity to upper case letters. 


genie . 

This 

etter? 


As an example of how to use this tool, a 51 x 21 90 degree 

duct is generated. The planes for this grid are: 


St reamwi.se 


1 

3 

5 

15 

25 

30 

34 

36 

37 
51 


( 1 -value ) 


Planes 

-6T5D 

-0 . 25D 
0 deg . D 
30 deg . D 
6 0 deg . D 
77.5 deg . D 
90 deg.D 
0 . 25D 
0 . 4 D 
2 . 5 D 


GENIE is a menu-driven program. The user defines the prob- 
lem domain through interaction with the program. To define the 
boundaries for a particular geometry, the first item (defining 
the boundaries) is selected. After this, the second item 
(creating the grid patches) is chosen in order to create the 
patches that are needed to generate the grid. In order to plot, 
the user picks the plotting item in the menu. All the data that 
will be required by the program in order to run makes up the date* 
file. By default, this is stored xn fort. 20. In order to pre- 

vent this file from being overwritten the next time a user runs 
GENIE, the file can be renamed or copied to another file using 
the UNIX command to move (mv filel fine2) or copy (cp 
f i 1 e 2 ) . The data file can be edited by using the editors 

are available on the workstations . Once the user finishes 


filel 
that 
the 
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interaction wi 
the problem that 
The output file 


the plotting 


routine, a g, ^phical solution to 
was defined by the user is viewed on the screen, 
(solution) is stored by default in fort 
is a binary file. An example of 

A. 


data file 


1 n 
can be 


47 
f ou nd 


This 
in Part 


t hat 

used 

age 

the 

order 


To plot the solution for the grid mentioned above, the steps 
were followed are located in Part B. The plotting package 

m^X 3 - 5 ' P1 ° t3d 3 - 5 i* • -"■> -d 3-D plotting ^ 

from NAoA/flmea. The next step that was taken was stacking 
c-D solution into 3-D. A gridstack program was written in 
-o accomplish this function. Since the output file that 


was generated by running GENIE was a binary file, this file had 

P? c ° nver ' tec * to a formatted file. In order to do this 

t lot 3d 3 . 5 was accessed again. The binary file was then read bl 
inputting ' > 


re/bin/x=fort .47 

fil e S is taX ° f the COmrnand to conve rt a binary file to a formatted 


list xyz/f ormatted /x=ou t pu tf i le 

. , . T |] e °^ tput f ile is specified by the user. This output- file 
will be the input file to the gridstack program. The gridstack 
TDAr^TL 15 called gridstack. f. The program is written in FOR- 
1 RAN . The program serves as a template in the sense that a user 
an ruH a ifferent g \ id simply by changing the dimension size to 
reflect the sice of the grid. The grid is stacked in the 7- 
direction using spacing from the Y-direction at inlet plane 
i ‘ * ^his program is located in Part C. This program was used 

to stack. ol x 31, 51 x 41, 101 x 41, 101 x 61, and 101 x 81 

grids. The gridstack program is compiled by using the compile 
command available on the IRIS. This command is 


f 7 7 f i 1 e name 


(in this case , gridstack. f) 


For t h i - 


To execute the file, the execute command is used, 
particular program, the execute command is 

gridstack.f -o gridstack 

For more information on these two commands, the system admini- 
trator should be consulted. Using the gridstack routine, the 

input file in this case was called induct. fmt and the output file 
outduct . fmt . 


After 
Plot 3d 3.5 
cuting the 
was 


compiling and executing the gridstack 
was accessed. Since the output file created 
gridstack program was a multigrid, the read 


p r og ram , 
by e x e - 
comm Cl nd 


re/ mgr /for / outduct 
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The sequence 
The different 


<- commands that followed are located in Part 
plots that are generated can be found in Part E 


D. 
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CONCLUSION 


GENIE was used to do geometry modeling and grid generation 
in support of engineer ing analysis at MSFC. It is recommended 
that the newer version of GENIE be implemented on the Silicon 
Graphics 4 - D system. This system is not currently available, but 
its acquisition is highly re com mended . 
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PART A 
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PROMPTING 

/Ito 

MUMBER OF POINTS IN THE I DIRECTION = 51 

number of points in the j direction = 21 

ACTIVITY OPTION Is DEFINING BOUNDARIES 
’NDICES ARE 1 51 1 1 FOR BOUNDARY 1 

— DUMBER OF SEGMENTS* 9 

1 3 ARE THE INDICES OF SEGMENT 1 OF BOUNDARY 1 

SEGMENT TYPE 2 STRAIGHT LINE 
LINE OPTION 1 USING 2 END POINTS 

INPUT POINTS ARE -.5000 -1.8000 AND -.2500 

PACK OPTION=0 EVEN SPACING 

3 5 ARE THE INDICES OF SEGMENT 2 OF BOUNDARY 1 

SEGMENT TYPE 2 STRAIGHT LINE 
LINE OPTION 1 USING 2 END POINTS 

INPUT POINTS ARE -.2500 -1.8000 AND .0000 

PACK OPTION=0 EVEN SPACING 

5 15 ARE THE INDICES OF SEGMENT 3 OF BOUNDARY 1 

SEGMENT TYPE 4 CIRCLE 
CIRCLE OPTION 3 USING 3 POINTS 
FIRST END POINT IS .0000 -1.8000 

OTHER END POINT IS .9000 -1.5588 

THIRD POINT ON CIRCLE IS .4659 -1.7387 

PACK OPTION=0 EVEN SPACING 
NO, GOING COUNTERCLOCKWISE 

15 25 ARE THE INDICES OF SEGMENT 4 OF BOUNDARY 1 

SEGMENT TYPE 4 CIRCLE 
CIRCLE OPTION 3 USING 3 POINTS 
FIRST END POINT IS .9000 -1.5588 

OTHER END POINT IS 1.5588 -.9000 

THIRD POINT ON CIRCLE IS 1.2728 -1.2728 

_ ACK OPTION* 0 EVEN SPACING 
NO, GOING COUNTERCLOCKWISE 

25 30 ARE THE INDICES OF SEGMENT 5 OF BOUNDARY 1 

SEGMENT TYPE 4 CIRCLE 
CIRCLE OPTION 3 USING 3 POINTS 
FIRST END POINT IS 1.5588 -.9000 

OTHER END POINT IS 1.7573 -.3896 

THIRD POINT ON CIRCLE IS 1.6751 -.6597 

PACK OPTION=0 EVEN SPACING 
NO, GOING COUNTERCLOCKWISE 

30 34 ARE THE INDICES OF SEGMENT 6 OF BOUNDARY 1 

SEGMENT TYPE 4 CIRCLE 
CIRCLE OPTION 3 USING 3 POINTS 
FIRST END POINT IS 1.7573 -.3886 

OTHER END POINT IS 1.8000 .0000 

THIRD POINT ON CIRCLE IS 1.7884 -.2038 

PACK OPTION=0 EVEN SPACING 
NO, GOING COUNTERCLOCKWISE 

34 36 ARE THE INDICES OF SEGMENT 7 OF BOUNDARY 1 

SEGMENT TYPE 2 STRAIGHT LINE 
LINE OPTION 1 USING 2 END POINTS 

INPUT POINTS ARE 1.8000 .0000 AND 1.8000 

PACK OPTION=0 EVEN SPACING 

36 37 ARE THE INDICES OF SEGMENT 8 OF BOUNDARY 1 


-1.8000 


-1.8000 


.2500 
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.2500 AND 1.8000 

9 OF BOUNDARY 1 


.4000 AND 1.8000 


SEGMENT TYPE 2 STRA1 T LINE 
LINE OPTION 1 USING 2 END POINTS 

INPUT POINTS ARE 1.8000 

PACK OPTION=0 EVEN SPACING 

37 51 ARE THE INDICES OF SEGMENT 
SEGMENT TYPE 2 STRAIGHT LINE 
LINE OPTION 1 USING 2 END POINTS 

INPUT POINTS ARE 1.8000 

PACK OPTION 3 0 EVEN SPACING 

INDICES ARE 1 51 21 21 FOR BOUNDARY 2 

NUMBER OF SEGMENTS 3 9 

* 3 ARE THE INDICES OF SEGMENT 1 OF BOUNDARY 

SEGMENT TYPE 2 STRAIGHT LINE 
LINE OPTION 1 USING 2 END POINTS 
INPUT POINTS ARE -.5000 -2.8000 AND 

PACK OPTION=0 EVEN SPACING 

3 5 ARE THE INDICES OF SEGMENT 2 OF BOUNDARY 

SEGMENT TYPE 2 STRAIGHT LINE 
LINE OPTION 1 USING 2 END POINTS 
INPUT POINTS ARE -.2500 -2.8000 AND 

PACK OPTION=0 EVEN SPACING 

5 15 ARE THE INDICES OF SEGMENT 3 OF BOUNDARY 

SEGMENT TYPE 4 CIRCLE 
CIRCLE OPTION 3 USING 3 POINTS 
FIRST END POINT IS .0000 

OTHER END POINT IS 1.4000 

THIRD POINT ON CIRCLE IS 
PACK OPTION=0 EVEN SPACING 
NO, GOING COUNTERCLOCKWISE 

15 25 ARE THE INDICES OF SEGMENT 

SEGMENT TYPE 4 CIRCLE 
- CIRCLE OPTION 3 USING 3 POINTS 
FIRST END POINT IS 1.4000 

OTHER END POINT IS 2.4249 

THIRD POINT ON CIRCLE IS 
PACK OPTION 3 0 EVEN SPACING 
NO, GOING COUNTERCLOCKWISE 

25 30 ARE THE INDICES OF SEGMENT 

SEGMENT TYPE 4 CIRCLE 
CIRCLE OPTION 3 USING 3 POINTS 
FIRST END POINT IS 2.4249 

OTHER END POINT IS 2.7336 

THIRD POINT ON CIRCLE IS 
PACK OPTION=0 EVEN SPACING 
NO, GOING COUNTERCLOCKWISE 

30 34 ARE THE INDICES OF SEGMENT 

SEGMENT TYPE 4 CIRCLE 

CIRCLE OPTION 3 USING 3 POINTS 
FIRST END POINT IS 2.7336 

OTHER END POINT IS 2.8000 

THIRD POINT ON CIRCLE IS 
PACK OPTION=0 EVEN SPACING 
NO, GOING COUNTERCLOCKWISE 

34 36 ARE THE INDICES OF SEGMENT 

SEGMENT TYPE 2 STRAIGHT LINE 


-.2500 


.0000 


-2.8000 

-2.4249 

.7247 -2.7046 


4 OF BOUNDARY 2 


-2.4249 

-1.4000 

1.9799 -1.9799 


5 OF BOUNDARY 2 


-1.4000 

-.6060 

2.6052 -1.0262 


6 OF BOUNDARY 2 


-.6060 

.0000 

2.7820 


-.3170 


7 OF BOUNDARY 2 


.4000 


2.5000 


-2.8000 


-2.8000 


16 



.4000 


LINE OPTION 1 USING 2 END POINTS 

P0INTS ARE 2.8000 .0000 AND 

PACK OPTION=0 EVEN SPACING 

SErMPN^rn^f ™ DICES 0F SEGMENT 8 OF BOUNDARY 

- SEGMENT TYPE 2 STRAIGHT LINE 

LINE OPTION 1 USING 2 END POINTS 

INPUT POINTS ARE 2.8000 .2500 AND 

PACK OPTION=0 EVEN SPACING 

qprMPum 55L TH V NDICES 0F SEGMENT 9 OF BOUNDARY 
SEGMENT TYPE 2 STRAIGHT LINE 

LINE OPTION 1 USING 2 END POINTS 

d»pvV° INTS ARE 2.8000 .4000 AND 

PACK OPTION=0 EVEN SPACING 

ii!SJS5 S rt ARE 1 1 1 51 FOR BOUNDARY 3 

NUMBER OF SEGMENTS = 1 

SEGMENT TYPE 2 STRAIGHT LINE 

LINE OPTION 1 USING 2 END POINTS 

D»«H T rt P0INTS ARE -.5000 -1.8000 AND 

OPTION=7 PACKED ON BOTH ENDS USING HYPERBOLIC 
WITH SMALLEST INTERVAL^ .0025000 

INTERVAL IN SECOND SECTION = .0025000 

INDICES ARE 51 51 1 21 FOR BOUNDARY 4 

NUMBER OF SEGMENTS = 1 

SEGMENT TYPE 2 STRAIGHT LINE 

LINE OPTION 1 USING 2 END POINTS 
INPUT FMNTS ARE 1.8000 2.5000 AND 

WITH QmIt?pcZ TMmE KED 0N B0TH ENDS USING HYPERBOLIC 
WITH SMALLEST INTERVAL= .0025000 

INTERVAL IN SECOND SECTION = .0025000 

0PTI0N 2s CREATING GRID PATCHES 

*** 1 51 1 21 F0R patch 1 

ACTIVITY OPTION 3: CREATING AUTOPATCHES 
ACTIVITY OPTION 5: PLOTTING 
1 PLOTS IN PLOT SET 1 
OUTER BORDER TO BE DRAWN ON PLOT 1 

1 PATCHES TO BE DRAWN ON PLOT 1 

1 51 1 21 FOR PATCH 1 ON PLOT 1 

ACTIVITY OPTION 8: OUTPUTTING FILE(S) 

YES, GRID IS TO BE OUTPUT 
YES, INPUTS ARE TO BE OUTPUT 


2.8000 .2500 

2 

2.8000 .4000 

2 

2.8000 2.5000 


-.5000 -2.8000 

TANGENT STRETCHING 


2.8000 2.5000 

TANGENT STRETCHING 
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PROMPTING 

vito 

JXJISpS ? F P0INTS IN THE I DIRECTION = 101 
NUMBER OF POINTS IN THE J DIRECTION = fil 

^NnrXi?e Y 0PTI0N 1 1 DEFINING BOUNDARIES 
NUMBER^OF^S egments ■ X , 1 1 ™ B ~ 1 

segment^type ™f 1 mSaa^S m 1 ° F B0UNDARY 

tmdtw°d TI0N 1 USING 2 END POINTS 

SEGMENT™ ™2 I s?^ E 1h? F L ?Se MENT J0PB0UNDARY : 

LINE OPTION 1 USING 2 END POINTS 

INPUT POINTS ARE - 2500 1 pnnn , Hn 

PACK OPTION=0 EVEN SPACING 00 1,8000 

SEGMENT^TYPE '^RCLE 0F SEGMENT 3 0F BOUNDARY J 

CIRCLE OPTION 3 USING 3 POINTS 

FIRST END POINT IS .0000 -1 8000 

OTHER END POINT IS .9000 -1*5588 

THIRD POINT ON CIRCLE IS 4659-1 7187 

PACK OPTION=0 EVEN SPACING 9 1,7387 

NO, GOING COUNTERCLOCKWISE 

SEGMENT ^TYPE ™4 * CIRCLE ° P SEGMENT 4 0F BOUNDARY 1 
CIRCLE OPTION 3 USING 3 POINTS 
FIRST END POINT IS .9000 

OTHER END POINT IS 1.5588 

HIRD POINT ON CIRCLE IS 
PACK OPTION=0 EVEN SPACING 
NO, GOING COUNTERCLOCKWISE 

_ ARE THE INDICES OF SEGMENT 

SEGMENT TYPE 4 CIRCLE 

CIRCLE OPTION 3 USING 3 POINTS 

FIRST END POINT IS 1.5588 - 9000 

OTHER END POINT IS 1.7573 -‘iflog 

THIRD POINT ON CIRCLE IS 1.6751 - 6597 

PACK OPTIONS EVEN SPACING 

NO, GOING COUNTERCLOCKWISE 

segment’type X c?rcle 0F SEGMENT 6 0F E0UNDART 1 

CIRCLE OPTION 3 USING 3 POINTS 
FIRST END POINT IS 1.7573 - 1886 

OTHER END POINT IS 1.8000 *0000 

THIRD POINT ON CIRCLE IS 1.7884 - 2038 

PACK OPTIONS EVEN SPACING 
NO, GOING COUNTERCLOCKWISE 

67 71 ARE THE INDICES OF SEGMENT 

SEGMENT TYPE 2 STRAIGHT LINE 
LINE OPTION 1 USING 2 END POINTS 
INPUT POINTS ARE 1.8000 

PACK OPTION=0 EVEN SPACING 

71 73 ARE THE INDICES OF SEGMENT 


-.2500 


.0000 


-1.5588 

-.9000 

1.2728 -1.2728 


5 OF BOUNDARY 1 


7 OF BOUNDARY 1 


.0000 AND 1.8000 


8 OF BOUNDARY 1 


-1.8000 


-1.8000 


.2500 
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.2500 AND 1.8000 

9 OF BOUNDARY 1 


.4000 AND 1.8000 


SEGMENT TYPE 2 STRAIGHT LINE 

LINE OPTION 1 USING 2 END POINTS 
INPUT POINTS ARE 1.8000 

PACK OPTION*0 EVEN SPACING 

73 101 ARE THE INDICES OF SEGMENT 
2GMENT TYPE 2 STRAIGHT LINE 

LINE OPTION 1 USING 2 END POINTS 
INPUT POINTS ARE 1.8000 

PACK OPTION=0 EVEN SPACING 

„ 1101 81 81 FOR BOUNDARY 2 

NUMBER OF SEGMENTS* 9 

?pr»Lwn 5 mS TH 5 INDICES 0F SEGMENT 1 OF BOUNDARY 
SEGMENT TYPE 2 STRAIGHT LINE 

LINE OPTION 1 USING 2 END POINTS 

-.5000 -2.8000 AND 

PACK OPTION=0 EVEN SPACING 

<?ErMPm* 5 ^ TH o INDICES 0F SEGMENT 2 OF BOUNDARY 
SEGMENT TYPE 2 STRAIGHT LINE 

LINE OPTION 1 USING 2 END POINTS 

ARE -.2500 -2.8000 AND 

PACK OPTION=0 EVEN SPACING 

INDICES 0F SEGMENT 3 OF BOUNDARY 
SEGMENT TYPE 4 CIRCLE 

CIRCLE OPTION 3 USING 3 POINTS 
FIRST END POINT IS .0000 

OTHER END POINT IS 1.4000 

THIRD POINT ON CIRCLE IS 
PACK OPTION* 0 EVEN SPACING 
NO # GOING COUNTERCLOCKWISE 

29 49 ARE THE INDICES OF SEGMENT 

3GMENT TYPE 4 CIRCLE 
'"'CIRCLE OPTION 3 USING 3 POINTS 
FIRST END POINT IS 1.4000 -2 4249 

EN 0 POINT IS 2.4249 -l^OOO 

THIRD POINT ON CIRCLE IS 1.9799 —1 9799 

PACK OPTIONS EVEN SPACING 
NO, GOING COUNTERCLOCKWISE 

INDICES 0F SEGMENT 5 OF BOUNDARY 
SEGMENT TYPE 4 CIRCLE 

CIRCLE OPTION 3 USING 3 POINTS 
FIRST END POINT IS 2.4249 -1.4000 

OTHER END POINT IS 2.7336 -.6060 

THIRD POINT ON CIRCLE IS 2.6052 -1.0262 

PACK OPTION=0 EVEN SPACING 
NO, GOING COUNTERCLOCKWISE 

59 67 ARE THE INDICES OF SEGMENT 

SEGMENT TYPE 4 CIRCLE 
CIRCLE OPTION 3 USING 3 POINTS 
FIRST END POINT IS 2.7336 

OTHER END POINT IS 2.8000 

THIRD POINT ON CIRCLE IS 
PACK OPTION=0 EVEN SPACING 
NO, GOING COUNTERCLOCKWISE 

67 71 ARE THE INDICES OF SEGMENT 

SEGMENT TYPE 2 STRAIGHT LINE 


-.2500 


.0000 


-2.8000 

-2.4249 

.7247 -2.7046 


4 OF BOUNDARY 2 


6 OF BOUNDARY 2 


-.6060 

.0000 

7820 


-.3170 


7 OF BOUNDARY 2 


.4000 


2.5000 


-2.8000 


-2.8000 
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LINE OPTION 1 USING 2 END POINTS 

INPUT POINTS ARE 2.8000 .0000 AND 2.8000 

PACK OPTION=0 EVEN SPACING 

71 73 ARE THE INDICES OF SEGMENT 8 OF BOUNDARY 2 

EGMENT TYPE 2 STRAIGHT LINE 
LINE OPTION 1 USING 2 END POINTS 

INPUT POINTS ARE 2.8000 .2500 AND 2.8000 

PACK OPTION=0 EVEN SPACING 

73 101 ARE THE INDICES OF SEGMENT 9 OF BOUNDARY 2 
SEGMENT TYPE 2 STRAIGHT LINE 
LINE OPTION 1 USING 2 END POINTS 

INPUT POINTS ARE 2.8000 .4000 AND 2.8000 

PACK OPTION* 0 EVEN SPACING 

INDICES ARE 1 1 1 81 FOR BOUNDARY 3 

NUMBER OF SEGMENTS* 1 
SEGMENT TYPE 2 STRAIGHT LINE 
LINE OPTION 1 USING 2 END POINTS 

INPUT POINTS ARE -.5000 -1.8000 AND -.5000 

PACK OPTION* 7 PACKED ON BOTH ENDS USING HYPERBOLIC TANGENT 
WITH SMALLEST INTERVAL* .0025000 
SMALLEST INTERVAL IN SECOND SECTION * .0025000 

INDICES ARE 101 101 1 81 FOR BOUNDARY 4 

NUMBER OF SEGMENTS* 1 
SEGMENT TYPE 2 STRAIGHT LINE 
LINE OPTION 1 USING 2 END POINTS 

INPUT POINTS ARE 1.8000 2.5000 AND 2.8000 

PACK OPTION* 7 PACKED ON BOTH ENDS USING HYPERBOLIC TANGENT 
WITH SMALLEST INTERVAL* .0025000 
SMALLEST INTERVAL IN SECOND SECTION = .0025000 

CTIVITY OPTION 2: CREATING GRID PATCHES 
'“INDICES ARE 1 101 1 81 FOR PATCH 1 

ACTIVITY OPTION 3i CREATING AUTOPATCHES 
ACTIVITY OPTION 5: PLOTTING 
1 PLOTS IN PLOT SET 1 
OUTER BORDER TO BE DRAWN ON PLOT 1 

1 PATCHES TO BE DRAWN ON PLOT 1 
INDICES ARE 1 101 1 81 FOR PATCH 1 ON PLOT 1 

ACTIVITY OPTION 8: OUTPUTTING FILE(S) 

YES, GRID IS TO BE OUTPUT 
YES, INPUTS ARE TO BE OUTPUT 


.2500 


.4000 


2.5000 


.-2.8000 

STRETCHING 


2.5000 

STRETCHING 
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PART B 


21 



In order to plot on Plot3d_3.5 , the following commands are followed: 

1 ) Access Plot3d_3.5 by inputting plot3d_3.5 at the 
prompt. 

2) Read the file that will be the input file to 
Plot3d_3.5. The format to read in Plot3d_3.5 is: 

re/ for or bin or unfor/x=name of the input file. 

The file created by running GENIE is a binary file. 

This means that using the above format, the read 
command will be re/bin/x=fort.47. 

3) Plot3d_3.5 allows the user to set the background of 
the screen during plotting. The syntax for this is bg 
#, where #, by default represents the number 

representations of the different colors that are 
available in the software package. The different 
colors that are available can be found in the 
Plot3d_3.5 user manual. I used a white background 
and the command is : bg 1 . 

4) At the prompt, enter wa. The following sequence of 
commands and input follows: 

Command Input 

wa 1, grid 1 

Enter I start ([,end],[,inc] ): a 

Enter J start ([,end],[,inc] ): a 

Enter K start ([,end],[,inc]): a 

Line Hidden, Line, or Shaded-Surface: 

Enter Color: rec j 

Enter Line Type: 

Enter Line Thickness Factor: 

Enter Symbol Type: 

Enter Symbol Size Factor: 

By inputting a , the user is specifying to the package that ail of 
the data in the I, J, K direction should be used. After going through 
these sequence, the user can plot by using the pi command. This plot 
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will be a 2-D plot. The plot is checked in order to make sure that it 
conforms to right hand co-ordinate system.. For more information on 
Plot3d_3.5, the user manual should be consulted. 
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PART C 
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c 

c 

c 


c 

c 


THIS IS A GRIDSTACK PROGRAM THAT STACKS A 2D GRID TO 3D 

THE n?i?f,!! SI0N SIZE REp LECTS THE SIZE OF THE GRID 

™ENE HM 
GENERATED BY PLOT3D USING THE LIST COMMAND 

open ( 5 , form= ' formatted ' , f ile= ' grid3 . fmt ' ) 
open ( 6 / f or xn= ' formatted ' , f ile= ' outduct fmtM 
READ IN 2D FORMATTED, SINGLE BLOCK GWD } 


READ ( 5 , * ) IMAX , JMAX 

I " 1 ' IMAX )' J’l.JMAX), L-1,2) 

DO 100 J= 1, JMAX 
DO 100 1=1 , IMAX 

X3(I,JMAXP1-J,1) = X1(I,J,1) + 0.5 

100 CONTINUE 1 ' JMAXPW ' 2 ’ = X1(I ' J ' 2 > + 2 ’ 3 

c INITIALIZE THIRD DIMENSION MAXIMUM INDEX AND BLOCK SIZE 

KMAX — JMAX 

c ( J- INDEX ^ AT InLt"plANE T (I= 1 ) K ~ INDEX) USING SPACING FR0M Y-DIRECTION 

DO 200 K = lficMAX 
DO 200 J = 1, JMAX 
_ DO 200 I * 1 , IMAX 

X2(I,J,K,1) = X3(I,J,1) 

X2 ( I / J,K, 2 ) = X3(I, J,2) 

200 X2 ( I ' J / R / 3 ) “ X3 ( 1,K, 2) 

200 CONTINUE \ t , , 

NGRID = 1 

WRITE ( 6, *) NGRID 
WRITE ( 6 , * ) IMAX, JMAX, KMAX 
DO 300 L = 1,3 

300 ^^CONTINUE ^ (X2(I ' J ' K/L ) '1=1/ IMAX), J=l, JMAX), K=1,KMAX) 

STOP 

END 
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In order to get the plotting lor the multigrid that is created by 

running the gridstack program, the following commands and input are 
done by the user. 

Commands , nput 


wa 

f 

a 

a 

red 


Plot3d_3.5: 
wa 1, grid 1 

Enter I start ([,end],[,inc]): 

Enter J start ([,end],[,inc]): 

Enter K start ([,end],[,inc]): 

Line Hidden, Line, or Shaded-Surface: 

Enter Color: 

Enter Line Type: 

Enter Symbol Type : 

Enter Symbol Size Factor: 

This sequence is repeated but this time the I, J, K value will be 

different. Also, the color will be different. The input are : 

I: Input a 
J: Input f 
K: Input a 
Color: Input green 


The sequence is repeated again. The I, J, K, and the color values are : 
I: Input a 
J: Input a 
K: Input f 
Color: Input blue 

After the last sequence, the solution is plotted by using the pi 
command. The plot will appear on the screen. The a input represents 
all of the input file data in the specific direction. The f input 
represents the first data in the specific direction. The I input 
represents the last data. The different colors are used for checking 
purpose. We want to make sure that the I, J, K direction are in the 
proper direction. Also, when the plot appears on the screen, the plot 
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is checked to make sure that it abides by the right hand co-ordinate 
system. 

The picture of the plot is taken by entering the T command at 
the prompt. After about 2 minutes, the screen will start to fade. A 
message "Screen Image is saved in outduct-a" will appear. The 
outduct-a file is the image file. Plot3d_3.5, by default, creates this 
file anytime the user uses the T command. The user needs to quit 
Plot3d_3.5 in order to send the image file to the printer. The quit 
command for Plot3d_3.5 is : quit. 

The user can also enter two lines of comment that will appear 
on the printout of the image file. The command for this is : t. This 
command is entered at the UNIX prompt not the Plot3d_3.5 prompt. 
After this, the image file is sent to the printer by entering the ftek 
command. The syntax for this command is : 

ftek name of image file. 

The name of the image file is the file that was created by Plot3d 3 5 
when the user took a picture of the plot (outduct-a). 
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OPTIMIZATION OF FORTRAN PROGRAMS WITH APPLICATION TO 
COMPUTATIONAL FLUID DYNAMICS 


INTRODUCTION 

Analysis of three-dimensional flow was thought impossible 
until equations governing the flow of air, water and other fluids 
were first published by Claude Louis M.H. Navier in 1823 and 
generalized by Sir George C. Stokes in 1840 ( Navier -Stokes equa- 
tions). The mathematics of these equations is too complex to 
allow an analytical solution in all but the rn os t simple cases; 
thus, numerical solutions using the most powerful computers 
(supercomputers ) currently available yield the best approxima- 
tions to the solutions of the Navier-Stokes equations. Conse- 
quently, numerical analysis and supercomputers have become essen- 
tial tools in Computational Fluid Dynamics (CFD) for the analysis 
of complex three-dimensional flows, ouch as the flow of fuel 
through the Space Shuttle Main Engine or the modeling of weather 
conditions around the earth. 

Most computationally intense CFD code is written in FORTRAN, 
although most current university curricula for engineers and 
scientists either contain no course work in FORTRAN or contain 
only a single elementary course in the FORTRAN language. Agen- 
cies such as the National Aeronautical and Space Administration 
(NASA) and the Department of Defense (DOD) as well as industry- 
have many large CFD codes which are often run and maintained by 
engineers and scientists who are lacking the experience or exper- 
tise of working with a FORTRAN code of this magnitude. Also, 
most have had no or very little experience with a supercomputer . 
As most of these users are not trained to modify FORTRAN codes to 
run efficiently on advanced architecture computers, they often 
have little knowledge or concern for the burden that these large 
codes place on supercomputer resources. Both the Cray X-MP at 
the Alabama Su pe r com pu t e r Center and the Cray X-MP at NASA/Mar- 
shall Space Flight Center (MSFC) are examples of supercomputers 
which are used to full capacity most of the time, 

BACKGROUND 

Many of the large CFD codes in current use by government and 
industry contain large amounts of FORTRAN written before 1980. 
At that time FORTRAN code was typically developed in a non- 
structured fashion using m any language constructs which are both 
outdated and extremely inefficient for modern supercomputers with 
highly developed compilers. These compilers can overcome some of 
the inefficiency in the FORTRAN code but leave a great deal of 
room for improvement by users. On current supercomputers, some 
mini-supercomputers (e.g. Alliant, Convex, etc.) and even on high 
performance workstations (e.g. Stardent , Silicon Graphics, etc.), 
improvement of the FORTRAN code's efficiency can dramatically 
improve performance. Improved performance yields reduced 
demand on already saturated computer systems, reduced design time 
for hardware and reduced customer cost for computing resources. 
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tbree moat effective means to improve code effi 
rori-,t'ion 9h vectorlzatior '> optimization and m icrotaski ng . Vec • 
* Pr ° C '” ^ -Me* * -ngle instruction is made ? to 

performed , ° Peratlo " s instead of Just one. This. process i. 

time°' of the c*. C 7 fcn ^pU l-'Wr'f 

important features o f ' *" 'IX “ v’i ^ \Z 

from central memory need to be made and arithmetic operation-'" 
simplified^ also reducing the computational time of the code 

solution 96 303 ff problems requiring supercomputers for their 
-olution are often constrained by insufficient power of < wr t 

L h :: e , h orde r : f r snitude - Modern -p^«.put e p r 

used i’ n .J 0le ? la f t0 dellver an order of magnitude more power if 
desirable 0 k T’ ^ maki "° optimisation of the code e 

tlt r : although , because of the inexperience of most codi 
*-t - ' lt - L - rarely done. Microtasking refers to the ahillf, jt 

the supercomputer to s imu 1 taneoTI^J 17 : 4. th bl> of 

on different centr-1 ™ b ' neous 1 V execute segments of a program 
e r - nt cent 1 o 1 processing units (CPUs) Mirrntv-i ; _ 

produces a code which allows the proiram -tr 1 ~ d 190 

one CPU i *• 1 . n . " e Progr am to be run even when only 

CPUS b<!COm ' -«il*ble the 

Mi Cro ta-Hnn H . ability to make use of the free CPU(s) 

the 't. ? d ° eS not r e< ^uce computational time, but can r-duc" 

cod r ihis ^ 

r t;m p ti„„ to m ak e optima! use^ ^^nt^n ^ cL “ h ° f 

relatively new feature and little pffnrf l ^ , 1 ^ 

CFD users on its use. U effort has oeen made to train 
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METHODS 


The objective of the work was to develop written guideline- 

and " ci r tiStS t0 and improve Te eTT 

^ o large scde CFD codes in current use or codes which 
under development for future use The wr ittpr n * u 7- 
to accompl i sh three goal,: -.ntten guideline, were 

1) Reduce the number of man-hour, required by the user who 
1 = lesponsible for improving the efficiency of- a code; 

2) Minimise the amount of control processing unit (CPU) 
time the code requires to run; 


3) Optimize the use of multiple CPUs available on 
supercomputer (parallel processing). 


the 


of Jftwarr m toolr' ne -!rren(I' e 0b:,eCtiv,, identification 

cuidel i npc a. H . r enlly in existence and development of 

FORTRAN code°f thS evaluation and restructuring of the 

current 1 -t-h ° r , S USer ' Cray Research Incorporated (CRI) is 

Suites V th ! t° n V m ^ ufactur ^ of supercomputers in the United 

run’rn’- ^‘ erefore ’ ail very large CFD codes are developed to he 

i " The T y C f mpUter * ThUS ’ the fir5t step toward accompl ish- 

fool- avail-ble W f f nCorporate the software development 

project' jt r int ° tHe Written sidelines of this 

MARK t - a T r f nt ^ ools available from CRI (FLOWTRACE, LOOP- 

MARK and Autotasking) are not readily used by scientist- and 

engineers since they require study and experience beyond the 

daiy scope of work of the typical user. The written guidelines 

+-|T U 6 TT K UC T° nS ° n h ° W 3nd when these tools can be used to 
the greatest benefit. ' w 


Before attempting code optimization 
ident i f y “the most: t; 
ce r .t rat e the efforts 
TRACE, which generate* 

1 ""rput trom 

by the program, the 

.. . . - r,e 5 t: be percentage of total 

tirn_- spent in each subroutine, the mi <mk*c - r+.e +- ; i_ . 

s u b r ou - 


t was important to 
.lfne-consuming portions of code and then con- 
on these. CRI offers a tool called FLOW - 

TpArc * . . , run-time statistics. Output from FLOW- 

TRA^E indicates the total computing time used 
amount of time spent in each subroutii 

ubroutine, the number of times each 


; ne ^ _ called , and the average time per call. It also produces 
o Railing tree which identifies the calls made at each level of 
:: e Program. FL0WTPACE w « used to determine which sections of 

. C ° e ~ lOU d investigated for possible improvement so that 
1 1 m © w o u I d n o t be spent 

little impact on overall 


improving a 
per forma nee . 


section of code which 


ha d 


Another tool available 
called LOOPMARK. 

FORTRAN code and 
vectorized by the 
each loop that is not 
which loops could be 


on the CRI line of supercomputers 


tion by the compiler. 


, ^ ^ _ w . -■ IS 

LOOPMARK examines each inner DO loop in the 
comments as to whether it is automatically 
compiler; it also prints an explanation for 
vectorized. LOOPMARK was used to determine 
restructured to allow automatic vectoriza- 
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Autotasking is a tool that is curr .:ly only available 
through use of the cf77 compiler. It first vectorizes the inner 
loop, then breaks up the outer loop, determining which portions 
Con be run in parallel, which variables c>re private and which are 
shared, and whether the loop performs enough work to warrant the 
overhead incurred by using more than one processor . The addi- 
tional overhead associated with Autotasking is incurred even if 
only one CPU is utilized; therefore invoking Autotasking when all 
CPUs are in use can actually result in an increase in CPU time 
with no decrease in wall-clock time. At least 85-70% of the code 
must be able to be run parallel for Autotasking to be beneficial 
j. n any case. Because the program statements ere no longer 
cuted in deterministic order, Autotasking car. also produce 
suits different from that of the original code. The user 
determine whether these differences are significant. 

OUTCOME OF TASKS 


exe - 
re- 
must 


Tm€? ou t come o f this work is 3 set oh written guidelines 
which will guide the code user through the efficiency evaluation 
of any CFD FORTRAN code. Thorough understanding of the mathemat- 
ics and physics in the code will not be necessary to carry out an 
analysis of the efficiency cf the code. Once the evaluation is 
^ ’ 1 r- 5 more guidelines are available to lead the engineer or 
scientist through an evaluation and modification of the code’s 
level of vectorization and optimization and the use of autotask- 
ing. This "tool box" of improvements provides easy to understand 
and easy to implement modifications to CFD codes which engineers 
and scientists can use to improve a code without having to know 
all the details of efficient FORTRAN programming. This set of 
guidelines will aid any engineering group whose main function is 
'"°de usage. New candidate codes can be immediately evaluated to 
determine their efficiency with respect to other codes in current 
use. Codes selected for use by the engineering group car, then be 


taken through an evaluation/modification phase 
demand on computing resources. 


resulti ng 
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GUIDELINES FOR EFFICIENCY EVA AT I ON 


Two tools, LOOPMARK and FLOWTRACE , are available from Cray Re- 
search to assist you in evaluating the efficiency of your FORTRAN 
code. LOOPMARK is run at compile time; therefore no additional 
overhead is incurred by using if. FLOWTRACE c ecu f os c.t i un 
time, and thus requires more CPU time than running the code 
without it; however, it provides valuable information that ran 
save you a great deal of time when trying to imprc.e the code’s 


ef f i ciency . 
or together 


LOOPMARK 


L OWTRACE 


eon b e u s e cJ e ithe r s e p a r c t e 1 v 


1* Compile the code using LOOPMARK . 
ing file with each DO loop mart 
toriced. Except f or outer loop 
iced, LOOPMARK lists the reasons 
i r ed . To compile y ou r 


LOOPMARK produces a list 
- d a e 1 t h e r scalar or \ e c • 
-■ , which cannot 


‘i at 


1 O O p 


c f t 


ode using LOOPMARK 
v ms g s file name . f 


oe .e c t o r ■ 
not vector - 
c i t her 


or 


. e na rue 


o r 

c f 7 7 ~ e rn f i 1 e n a m c . f 

The following is an example of output from LOOPMARK: 


n o 


ORIGINAL PAGE ?S 
OF POOR QUALITY 



175 

176 

177 

178 


175. S 

176. : v < 

177. v 

178. : V > 


DO 50 K-l , K£ND(NT) 
DO 50 I-l.IMXI(NT) 
QW(I,K)=HF 
50 CONTINUE 


VECTORIZATION 


INFORMATION 


• Loon !f &r ^ ne at line 124 was vectorized 

• Loon at Hne 128 was vectorized 

• Loon at line 133 was vectorized 

the P in Ur at line 139 was not vectorized 

, i™ lo °P contains input/output operations 

th P . S arting at line 144 was not vectorized 
„ f the loop contains input/output operation 
Loop starting at line 151 was not vectorized 
he loop contains input/output operations 
Loop s Ur ting at line 160 was not vectorized 
lo fP contains input/output operations 
th P , s arting line 164 was not vectorized 
tnn °? P c ° nta * ns input/output operations 
th P . s artlng at line 169 was not vectorized 
the loop contains input/output operations 
Loop starting at line 176 was vectorized 

th P . s Artin * at line 233 was not vectorized 
oop contains a scalar store on NTOTAL 


because 

because 

because 

because 

because 

because 


because 
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Run the code using FLQWTRACE to determine in which subrou- 

FLOUTRA^" m °^ eff ° rt ChOUid bC conc * ntr.ted. Output from 
FLOUTRACE indicates the total computing time used by the 

program broker, down into the time spent in esc h subroutine 

^ nc ' ^ e percentage: of totc*l time “pent in each. It 

indicates the number of times each subroutine is called fd' 

the average time per call, and produces soiling tree which 

identifies the calling routine for each. When attempting to 
improve the code's efficiency, you should first examine 
,'JtT subroutines in which the largest percentage of the 
,r’; al tlrne AS --pent. You should also examine ar.r subroutine 
-hat is called a large number of times. Although the per- 
centage of total time spent in it may be small, improving n 
subroutine that i, called many times can decrease the total 
b,- J time required to run the code. FLOUTRACE can be invoked 
using either the eft, cft77 or cf77 FORTRAN compiler bv 
using the following compile commands: 


f t 7 7 - e f f ilen a me . f 


cf 7 7 -F 


The foliowii 


example of output from FLQWTRACE : 
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1 flowtrac 

E — Alphabetized 

summary 

0 Routine 

Time executing 

Called 

24 BTRI 

6.632 

( 14.40*) 

510 

©00062361a 

10 CONTRL 

0.438 

( 0.92%) 

5 

43 CONVRG 

0.245 

( 0.52*) 

5 

39 CORREC 

0.419 

( 0.88*) 

5 

20 FLUXR 

3.300 

( 6.95*) 

37360 

©00047205a 

8 CETAJA 

> 

( 0.00%) 

13 

©00111741a 

14 CETOLD 

> 

( 0.00*) 

20 

•00111631a 

7 GETQ 

0.002 

( 0.01*) 

70 

©00111565a 

33 GETTMP 

> 

( 0.00*) 

5 

9 GETXYZ 

0.004 

( 0.01*) 

40 

©00112004a 

46 INEXS 

0.005 

( 0.01*) 

1 

22 LHSSI 

3.190 

( 6.72*) 

5 

35 LHSSO 

1.755 

( 3.70*) 

5 

18 MUKN 

0.288 

( 0.61*) 

10 

©00047027a 

19 MUTRSI 

3.439 

( 7.25*) 

5 

34 MUTRSO 

1.070 

( 2.25*) 

5 

2 OPNCLO 

0.009 

( 0.02*) 

1 

49 OUTPUT 

0.723 

( 1.52*) 

1 

8 PUTAJA 

> 

( 0.00*) 

3 

13 PUTOLD 

> 

( 0.00*) 

10 

5 PUTQ 

0.001 

( 0.00*) 

38 

©00111216a 

30 PUTTMP 

> 

( 0.00*) 

5 

4 PUTXYZ 

> 

( 0.00*) 

3 

3 REDWRT 

16.794 

( 35.39*) 

2 

17 RHSSI 

1.082 

( 2.28*) 

5 

26 RHSSO 

0.840 

( 1.77*) 

5 

1 ROTOR 

0.003 

( 0.01*) 

1 

23 SMATRX 

5.092 

( 10.73*) 

37360 

•00061121a 

31 SRINT 

0.016 

( 0.03*) 

5 

21 VFLUX 

0.921 

( 1.94*) 

26090 

•00046571a 

25 VMAT 

0.974 

( 2.05*) 

26090 

•00060621a 

* • * TOTAL 

47.446 

127683 Total call 


1 FLOWTRACE — Calling tree 
1 ROTOR 00025203a 


2 

OPNCLO 

00111164a 

3 

REDWRT 

00106441a 

4 

PUTXYZ 

00111435a 

5 

PUTQ 

00111216a 

6 

PUTAJA 

00111372a 

7 

GETQ 

00111565a 

8 

GETAJA 

00111741a 

9 

GETXYZ 

00112004a 

10 

CONTRL 

00027747a 

11 

PUTQ 

00111218a 

12 

GETQ 

00111565a 

13 

PUTOLD 

00111262a 

14 

GETOLD 

00111631a 


Average T 


013 Called by 

LHSSI 

LHSSO 



315 

195 


088 ©00027747a 

Called by 

ROTOR 


049 ©00031204a 

Called by 

CONTRL 


084 ©00031526a 

Called by 

CONTRL 


> Called by 

RHSSI 

RHSSO 



23395 

13965 


> Called by 

CONTRL 

REDWRT 



10 

3 


> Called by 

CONTRL 

CONVRG 



10 

10 


> Called by 

CONTRL 

CONVRG 

CORREC 


20 

10 

30 

> ©00111675a 

Called by 

RHSSO 


> Called by 

CONTRL 

CORREC 

INEXS 


10 

25 

2 


0.005 ©00026422a Called by ROTOR 
0.638 ©00047747a Called by CONTRL 
0.351 ©00054126a Called by CONTRL 
0.029 Called by RHSSI RHSSO 

5 5 

0.688 ©00063470a Called by RHSSI 
0.214 ©00077761a Called by RHSSO 
0.009 ©00111164a Called by ROTOR 
0.723 ©00026072a Called by ROTOR 

> ©00111372a Called by REDWRT 

> ©00111262a Called by CONTRL 

> Called by CONTRL CORREC REDWRT 

10 25 3 

> ©00111326a Called by RHSSO 

> ©00111435a Called by REDWRT 
8.397 ©00106441a Called by ROTOR 
0.216 ©00035550a Called by CONTRL 
0.168 ©00040363a Called by CONTRL 
0.003 ©00025203a Called by 

> Called by LHSSI LHSSO 

23395 13965 

0.003 ©00045133a Called by RHSSO 

> Called by RHSSI RHSSO 

21210 4880 

> Called by LHSSI LHSSO 

21210 4880 



IS 

GETAJA 

00111741a 

16 

GETXYZ 

00112004a 

17 

RHSSI 

00035550a 

18 

MUKN 

00047027a 

19 

MUTRSI 00063470a 

20 

FLUXR 

00047205a 

21 

VFLUX 

00046571a 

22 

LHSSI 

00047747a 

23 

SMATRX 

00061121a 

24 

BTRI 

00062361a 

25 

VMAT 

00060621a 

26 

RHSS0 

00040363a 

27 

MUKN 

00047027a 

28 

FLUXR 

00047205a 

29 

VFLUX 

00046571a 

30 

PUTTMP 

00111328a 

31 

SRINT 

00045133a 

32 

GETQ 

00111565a 

33 

GETTMP 

00111675a 

34 

MUTRSO 

00077761a 

35 

LHSS0 

00054126a 

36 

SMATRX 

00061121a 

37 

BTRI 

00062361a 

38 

VMAT 

00060621a 

39 

C0RREC 

00031526a 

40 

PUTQ 

00111216a 

41 

GETQ 

00111565a 

42 

GETXYZ 

00112004a 

43 

CONVRG 

00031204a 

44 

GETQ 

00111565a 

45 

GETOLD 

00111631a 

46 

INEXS 00026422a 

47 

GETQ 

00111565a 

48 

GETXYZ 

00112004a 

49 

OUTPUT 00026072a 


LOOPMARK and FLOWTRACE can be invoked simultaneously by 
using the commands: 

eft -V msgs filename. f 
eft 77 -e f m filename. f 
cf77 -F -em filename. f 


NOTE: 


It is advisable to run the code after any major 
changes are made to insure that the output is the same as 
that of the original code. If several changes are made 
before the code is run, it can become extremely diffi- 
cult to determine which change caused the discrepancy in 
t he ou tput . 
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GUIDELINES FOR VECTORIZA1 N 


A loop containing 

a no t he r 

loop 

will 

n e v e r 

vector 

i ze . Be - 

c c u s a t he inner rn o s 

t loops c 

a n b e 

v ,'-rt 

c r i r : 

d , the 

s f f L c i e n c > 

of the code c a n 

be greatl 

V i m p 

nr-' e d 

b y e 

n suring 

that tine 

1 1 'i n e- ; rn ost loop i s 

e x ecuted 

f o r t 

he 1 

rgest 

n urn be r 

of iters 

tions. A section 

of code 

c o n t a 

i n I ng 

nes t 

c d loops. 

with nc 


s t o ternentLi between . LSI 


Example: DO 10 I = 1, L 

DO 20 J = 1, M 

DO CO Y = 1 , K 1 

c a n I ways be r e o r d c* r e d w I Lb t ■> ^ . > n e ■ p e c t c* d s f f o c t . 
other statements occur within the ‘.ester! loops 

Example: DO 1 0 I - 1 , L 

s t a t erne nt s 
DO CO -! - 1, M 

s t a t c: rn c n t c 
DO 3 0 K = 1, N 

the loops 
rn o v e a n y 
happen! n g 

2 . An outer loop which is executed fer four or fewer iter .. ~ 
tions , e v e n I f v ectorized by the compiler , ca n be 
"unrolled 11 , which may significantly decrease the amount of 
CPU time used. (This is somewhat machine dependent.) The 
term "unrolling 11 means to repeat that section of the code 
for each Iteration, putting in constants for the array 
subscripts Instead of variable names. 

Example: DO 10 J = 1 , 3 

DO 10 K = 1, 1000 

10 A ( J , K ) - expression 

c a n b e u n r o 1 is c J t o produce 

DO 10 Y = 1, 1000 

A ( 1 , K ) - e x p r e ssio r- 

10 A ( 2 , K ) - expression 

0 . A loo p w i t hi o n 1 t w o iteration s will n c t b c v c c t o r i z o d ! :< : . 
the compile r a n d z h o u 1 d a 1 w a y s b c unrcl 1 c d . 

4. An inner loop with few iterations may be marked b> LCOPMAPP 
as a M Vs" (short vector) loop. In this case, no benefit can 
be produced by unrolling the loop. 

5 . An y loop that contains a READ or WRITE s tatement will n c t 
v e ctori z e . A WRITE s t a t e m e n t w i t h 1 r. a r ( IP s t a t c rn c n t , such 


can still be reordered , but you must be carr^vl 
■ t a t e n i e nt s in b e t w e e n t hi a t m a y c h a n go w In o t 
i n the code. 


d 3 



Ih DEBUG THEN 

WRITE (B, 1000) 

t o 1 1 f - 1 ' h 1 1 " 1 '' t T C ! ° ' 1 ^ “ 1 1 ° n even if the value of DEBUG i - , e t 

r ? ?-»■: -n- 

it), v.eCr-,i-; : ep T ate , °° p bef ° re th ° 100,1 «*.. 

' 1 veci -^ ( uaUoh can be achieved. 

A loop that contain*: a rail •* - ... . 

user-defined f,,,--- " 1 "' J " J - u br ou 1 i ne or eternal 

f jiil i. ion i>Ji 1 1 not vector i -;c- Mn -• ■ 4 . • . 

functions (SORT, SIM etc 1 - r - " ir, - lc 

1 1 * 1 V L v. . J aTC VC'C t Cl' r ) *•* A r r- 1 ‘ ■ _ 

.«:=ri-^- u ,c th ~ * r * f =°": 

vectot ization are: " IC '"" 

c“l* fl,ncti ° r ' * wt «i Of -ur, external function 

' ( NOTE . Be sure to look at the cal linn h-r C 

S“ „ l ! ,. Fl :T RACE determine uhich other Lbrjl* 


d . 


. ^ u-lr er sui 

f u^-- i! * funCtiDn if you decide to eliminate 
function ci 1 together . ) 


+ ht 


Move the entire subroutine code into th* loo- ( NOTF 

rp. r l Ur ° t0 r 001 ?* the caliin, tree „r ^ ' Vv ! 

■RpCE to determine which other subroutines .-,,11 7-it- 

subroutine if you decide to eliminate the i r,r 

Move the loop into the subroutine and moke , call r 0 ■ - 
1 om The original subroutine. 


This requires a more thorough understanding r 

d d C . ar l only be dor,e if The called s.ubrouti 
called from any other part of the program. 

Move the call outside the main loop and into 
loop of its own. This can only be done if th< 


' V t t z. 


1 <c t 


. tioni 


ire i ii e t : 


-■e pc. rate 

t Q i lowing 


Th,e_ called subroutine does not make 
variables referenced in the loop. 


. 2 n m c n t 


There are no STOP or alternate RETURN - taterncr.t- 
the subroutine. - 


The subroutine arguments are not arrays 
Any subroutine called from within this 
meets these same conditions. 


.subrou t i nc 


loop containing ® backward GO TO will not .ectorio- r ,i 

Ieoh' r .h t ^ ,CtUre h Structuring requires analysis ’of ‘ °hc 
1 . - ,no I* -cmecirnes impossible. When it is possible it 

requires that the loop be rewritten with a forward GO TO or 
- jump cut of the loop when a desired value is , cached (see 
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# 9 ) . 


2. A loop containing on conig nod GO T 0 w i 1 1 not vectorire nd 

cannot be restructured without rewr it ifig tc _ 1 irnir,; te the 

A S S I G N statement, which n ■ a y n o t b c p o c c i b I e i *,< i t h o u t a t h o r - 
o u g hi u n derctanding of t hi e code. 

9- A loop containing a jump out of the loop based on on IF 

ctatemen t w ill not / e ctorizc . T \ \ i c c a n c o rn et 1 me s b e r w f ■ i t. - 

ten so that the rest of the loop is only executed if the 
opposite of the IF condition is true. (NOTE: "“hi-., .lot: not 

a 1 w a y s s p eed up e.^ec u t i c r i . ) 

Example: DO 10 I - 1 , M 

stateme n t c 

IF (A . GE. B) GO TO 10 
e t a t c rn e n t s 
10 CONTINUE 

can be rewritten as 

DO 10 I = 1, N 
s tate rri c r, 1 s 
IF (A . LT . E ) THEN 
stateinen t s 
END IF 
10 CONTINUE 

In cases where this simple technique will not wort, there is 
no easy way to restructure the loop without . complete 
understanding of the code and a thorough t now ledge of FOR- 
TRAN. Those interested should consult the discussion jf 
stripmining in A Guidebook to FORTRAN on Supercomputers by 
John M. Levesque and Joel W. Wiiliarnson. 

10. A loop containing two IF statements with mutually exclusive 
c o n d i 1 1 o n s t h a t h a v e c c alar v ari a b 1 es be i n a update d d e p e r d - 
i ng on which condition is true will be flagged by LOQPM-PF 
as non - vector izable because ’’scalar- -.lues arc updated v 
than once' 1 . The loop will vector ice If It is rewritten 
using an I F ~ T HE N - E L S E t r u ctu r e . 

Example: DO 10 I - 1, N 

s t a t e m ents 
IF (A . NE . B) THEN 
statement set 1 
END IF 

IF (A .EG. B) THEN 
at a t e m e n t s e t 2 
END IF 
10 CONTINUE 

c a n b e rewr i 1 1 e n a s 
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1 1 . 


DO 1 0 I = l , N 

IF (A . NE. B) THEN 
statement set 1 
ELSE 

s t a t e n i e n t : c t D 

END IF 
10 CONTINUE 


A 

loop con 

tain 

, i Pi 

vec 

t o r i o e . 

Th 

i s 

on 

the right ha 

nd 

left hand a 

ide , 

a 

i s 

the resu 

It o 

r 

t h e 

a r ray 

w i t h 

C2 


: " c c lj r s i o n ~> r, n r r a y o r 

L-l S 1 Pi y cj : Cci 1 -.'i P O 

:jup t. 1 or, be fore I r ,g 
nd 5 - sig n .i og a n .. r r -jy c 1 une n t 

-■ •' 1 c u 1 a 1 i o ; , J n v olving a n at!' 


j f .a r 


expression . 


'■^ubijcr i pt clcci" Ffii e n ted b y 


- c, -* ; i -s. r ^ ill not 

ai~r ay conol ant 
t c n t In e 
value that 
element of 


. o n s t a r, t 


E x 3 rn pie: DO 1 0 C I = ? N 

100 A(I) - A ( I ' - 1) + e * pi-i 


In this case, LOOPMARK will flag the loop with « .. r .,- u .- r , 
was found on A” This car, be resolved by c taring fbpv' 
* in a temporary army. 


of 


do i o o i - i 1 rj 

100 TEMP ( I) A ( I 

DO 000 1 = 2, N 

2 0 0 A ( I ) = TEMP ( I - 1 ) + e p ,, £ ,, c . r , 


Although this permits 
b r f e a *3 i b 1 e be c a u s e c f 


/ector i ration of the loop. It may not 
t hi e a d d i 1 1 o n a 1 rn c rn o r y v ~ p q l i i r e rn c r, t 
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GUIDELINES FOR OPTIMIZAT >1 


If the compiler can recog nice in.ai i ant code within th 


loop, it can pre -compute it 
result i n a r e gist c r , w h i c I 


:e " 


; r o the- loop a nd s t o r e 


.an greatly reduce CPU t i m e 
P 1 a c i n g the i n v or i a n f p o r t i o n o f code L n side p a r e n t h c s c 

a *e e 1 1 . t s t h e c orn pi 1c r 1 n t h e t • c- cog n ; tion p r o c e s a . 


E >: a rn pie : 


DO 1 0 I 


1, N 


10 A ( 

I 1 - X » 

3 ( I ; * V 




is optima 1 1 y w r itte n 

o S' 





DO 10 

1 = 1, N 





10 A ( 

T ' - < v i 

V) ^ EC I) 




(NOTE. With the c,p 

mo r e o f f icient t ha n 

abi 1 I t i os 

of current ~ 

orn p i 1 ~- 

rs , t 

h i s 

TEMP 

- X +• Y 





DO 1 0 

I a l , ri 





10 A ( 

I .) - B( I ) 

♦ TEX-' 

) 



S u b c x p ressions cor n rn on to two or m ore e x 

pres a- i 

o n s 

c :> i i 

computed once and 

stored i n 

registers if 

the 

c om p i 

ler 

a b 1 c to r ecog nice t he rn . T he 

compiler 1 s effort 

can 

b e 

ha need by placing t h e c ube x p r 

e- ssion in par 

enth e *s 

es 


E> amp 1 e : 






A ( I ) 

= B(i) + 

( C ( I ) * D(I)) 




E( I ) 

= F ( I ) + 

G(I ) 




H(I) 

- P ( I ) ♦ 

(C(I) * D(I)) 

/ QO 

> 


A g a i n , t h i s i . ' rn o r e 

efficient 

than placing 

t h e ■ 

o m i n o n 

1 u 

pres i ion i ^ j> t c r n p o r 

1 ' ■ ■/ a r i z> 

bit a nd u s i ng 

the t 

empeo- 

a r c 


] f t w o 3 hi o r f loops contai n c e m m o n s ube pro s c i o n c 
1 o o p l c a n b e combined w i t h o u t c I . a r , g x n g t In e rn e a n I n g 
code, do so. This allows better uti lizct io n 

c o rn \ > x 1 c r ‘ r op t i rn 1 z a t icn c a p a h i 1 1 1 i e s . 


of 


t h 

4- k 


Mul tiplication operation^ are less -omputat ionn 1 iy s-pensiv 
t ! i a n a ; “ e d i v i s 1 o i> o p c i ' c tion s . C h a n g e d L v i s> i o n b y a c o n s t a. i n 
to multiplication whenever possible without losing accuracy. 


E x a m p 1 e 


A(I) - B( I ) / 


A ( I ) - B(I) 


4 7 


becomes 



Convert floating point exponentiation integer exponentia- 

tion whenever possible . 

Examp 1 e : 

y * * 4,0 

becomes 


X * * 4 


48 



GUIDELINES FOR MICROTASK i 


Micro 
the code 
d irectiver 
sired. Us 
that the 
u nders ta nd 
tlon ■?. t’d 
different 
■: . - d ; n t 


tasking can be accompli s h e d auto mu t i c c 1 1 y t h t " o u g h o u t 

i j s i n g a u to tc. s king or s e 1 e cti v e 1 y b y inc 1 u d i r i g c ompi 1 e r 
i n t he s pe c i f i c s ection s w h e: r e rr# i c r o t a o k i r.g is cJ e - 


e or compiler directives to force rn i c r ct a r- k i ng 
user be thoroughly familiar i^ith the cede, 
the scope of the variables affected by the code 


r equi r e s 
He mus t 


! j ride r s t a n d the effect s o f c I n i j 1 
CPUs. CPI provides fools such us 
bo d e t e r n i s n < t i c , n of scope. T I 


: a ' 'i o o u _ 
SPY and 


,-rf , 


e i n t c r o e t c d 


compiler d Ire c t i v e s t o enforce microt j s k I n y s houlc! r : for J ; o f h c 
CRI User 1 o Guide to cf7 n . 


Autotasking may be i n v o k e d u s i n g t h e comma n c 

cf 7 7 -ZP 


■ jn d may v 1-0 be used in combination with LOOPMAPK and/or FLOW 
T P A C E w i t h t h. e comm a 1 1 d s 


c f 7 7 -err. -Zp (LOOPMAPK) 

cf 7 7 -F -Zp (TLOWTRACE) 


cf 7 7 -em ~F -Zp 


( b c- 1 hi ) 


Autotas 
required to 
available , 
total number 
u s e f u 1 w h e n 


b 1 n Q does not always reduce the 
run the code because it depends on 
r i d c a n o nig. reduce the time by 1 / n 
of CPUs, at best. Autotasking is 
running the code at off-peak times. 


wa 1 1 - c 1 oc k 
the number of 
, w In e r e n i s 
more iikelv t 


t i rn e 
CPUs 
the 
be 


Be sure to examine all output after running the cod 
autotasking, as statements may be executed in a different 
x n d rn a y produce different results 


i th 
d e r 
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RESULTS 


These methods and techniques were applied to a 
large CFD code called ROTOR, which was obtained f 
shall Space Flight Center. This version of ROTOR 


verrion of a 

roni NASA /Mar " 


routine, including the main program. Of 


w h i < 


these 




were called relatively few times and in which G.OlA 
o f t h e total C P U 1 1 rn e w a s s pent. Of t he r c- a ; I n i n g 1 3 

tines, four were completely vector iced , including the 
w 1 1 j. c. [ * i the- 1 a r g e s t p e r c e n t 3. c; c s c f t i rn e w e r e r p e n t , a n d 1 i 
in put /out put routines. Therefore, - f forts to tr.hcnc.e 


‘re w e i ' 






jbre-u 
wo i r i 
i jere 


i rna r i 1 y 1 j_ rn i t e d 


of 


. n e : 


Vector i zat ion 


e n \ 


i n c e m e n t 


rf forts in three of thee 
tines resulted in an increase in the amount of CPU 
largest increase was 8.988 units, with ?» mean increase 
units. The unit time deer e a. s e obs e r v e d in 1 1' t e o t h c r s 
ranged from 0.082 to 58.771, with a mean decrease 
units. When efforts in these subroutines were exhauste 
j.n .significant subroutines were also o x a n i i n e d a n d 

After all modifications had been made, the total time f 
tion j. f the program dropped from 1 058.703 to 981.299 
decrees 1 


the program dropped from 1058.703 to 961. 
of 97.415 units., or 9.2V 


e s 
1 1 m e 
of 

u b r c> 

of 

d , t 

ill O d 

c*r 
u r i i 


U u f OU ’ 

; the 
3.077 
uti nes 
11 .422 
he 13 
i f i ed . 
e ecu- 
tr- . a 


c- ami no tion of each subroutine revealed that no 
sate, on could be perf orrned . 


c o d c _> p t i m i 


compi led a ncl run w 1 1 \ 
for output parameter 
ode in the 


The code was 
e n a Ij 1 e d . T h e v a 1 u e 

o u t put of the original 

Tot DROMAX in t It e s i x t h decimal place. Inlet 
results differed from those obtained from the 
the code as follows 

Parame ter* Decimal pi 


the autotask ing f e a t u r e 
DPIMAX differed from the 


n i n t h decimal place 

T ^.4- - nd r- y ] f 

- 1 


a nd 


' i g s n.: 


RO I NL 
RUINL 
RVINL 
RWINL 
PS I NL 
PTINL 
ROEXT 
PSEXT 
PTEXT 
V e 1 o c i t 


t ha t 

4 : i o n 
:-r, of 


The 
c a n 1 


e r n i u s t d e t e. i “ rn i n e w h e t }"» c r t h c s c differ e n c e s 


c i g n i f i 


ORIGINAL lb 

OF POOR QUALITY 



APPENDI 
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I mplementat ion of Equilibrium Chemistry 


1- Introduction 



The 

k now 1 

edg 

e of 

c hem i ca 1 

equ i 1 i br 

i u m c orn poo i t 

I o n s 
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c hem i 
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nozz 1 o: , 

a ng i n os , s I' i 

Gck 

tubes , 

heat 

e x c h 

a nge r s 

i u 

n d p r 

ocess i n g 

equ i pine n t 

. 




Considerable numerical calculations a re necessary to obtain 
equilibrium compositions for complex chcmir 1 : y a t oms . This has 

resulted in a number of digital computer program 1 : to do the 

calculations. A computer program (CEC , 1) written at NASA Lewis 

Res es rch Cc ntc r in 19 61 - 19 6 2 has had a w 1 d e accepta nee . Howe ver , 
1 n ca leu 1 a 1 1 ng t he f a s t hem leal ; cacti o ns in large-seal e a nd 
mu 1 1 i d i rvi e n sic n a 1 fi r o b 1 z m r- , t b: is ge n c r a 1 c fa emis t r y code (CEC) i c 
so inefficient. There arc sever.:. 1 r equ i of n c qu 1 1 1 b r i. ■ im 

chemistry solver to be implemented in CFD codes. These- art the 
computational efficiency and the reliability. T b.e main objective 
of the present study is to implement a n d modi f y the e x i s ting 
equilibrium chemistry solver to meet these requirements . 

To improve the reliability, the equilibrium solver- with the 
free energy minimization p r occ-c! u 1 * : is s e 1 e c t ed . C h e m i c a 1 

equilibrium is usually described by cither of two equivalent 
formulations - equilibrium constants or minimization of free 
energy . In comparison with Gibb's free energy minimization 
method, the equilibrium constant method has several disadvantages 
which are mo r e b o o !- V n e p 1 r i g , in u m c r i c a 1 difficulties:- with use of 
chemical components , and more difficulty.- in testing for presence 
of some condensed species.. Furtherniorr , f he Gibbs free energy is 
most easily minimized. Fo:~ these reason: , the free- energy 
m i n irniz a lion f o r m u 1 -a t i o n i s u s c d . 

T he p r i rnary requ i r erne nt of a n e q u i i i b r i u m s o 1 v c- r i n a 1 a r g e 
f i n i t e - d i f f e r e n c e *: o d e is t h a. t i t b c f a st . A ty p i c a 1 r u n o f 
multidimensional t r a n * i e n t e pre y r c a t, t i ng f 1 o w s r e q u i r e s t h c 
equilibrium p r oble rn t o b e a o 1 v e d b o u t a m i 1 1 i o n t i m e s: . 
Typically , the speed of an algorithm in solving a single problem 
i s n t easured i n rn i lliseccndo . Esc h a >•: t r a m j 1 1 is e c o n d r e q u i r e d b y 
t h e sol v e r w ill t raivs la te into about 19 a. ddjtional min u t o s o f 
corn pu t e r CPU t i me . To impr ov e the e f f i c i r nc y 1 r i t c r ms of c orn pu f ~ 
i ng time arid c o n v e r g e n c e , t hi e p r e s e r> t e q u i 1 i b r it: m c I - 1 e m i s t r y 
s o 1 v e r i ^ i 1 1 a d o p t : hi i e r a t ‘ c h i c a 1 algor i t h rn [ 2 ] w h * c h co rn Is i n c t h c- 

speed of Gauss ■ Siedcl method ?nd the convergence of Newton meth- 
od . 
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2. General Progr 


Features 


The present cqu i 1 i b r i urn code (ADAPT) has the f u nc t i o na 1 1 j 

similar features wit h o t h e r e x i s t i n g s q u 1 1 1 b : : urn c o 1 . o r c •" * c c p t 
t he adaptabi L i t y j i t h t h e m u It id imens io n a 1 C F D c o d e e . ADAPT 
eolation procedure using C. i b b c ’ tree c n e r g y n i s. n i rn ica h i o n 

t g c h niqu e p r o v ides useful opt i o n s sit h g r' a comp! c t c * cembu s t i o n 

equilibrium model or a full e q u i 1 i b i ' i u rn m g* del, in w h ( irh t h s f u 1 1 
: r r c.y o f s pec i cs c r i s t i ng u ndtr c hern i c a 1 eq u i 1 i b r i u m c ondit ions 
is obtained at each grid point with the given pressure nd 
e rit ha 1 py . 


In gene r a 1 , t hi e A D A P T code provi J g - 
ms t hemat Leal mod e 1 o f t be r arnj e t c orn be c t g ; 
be used to estimate the performance g. 


o rn p u ter -orient e ci 
n c c c: 1 c w h i c h c a n 
■ a m j e t d esig n a nd 


p a r a me t r i colly e v a 1 u a t e t he c f f c. c t s of c I a t ig o ■ s in the design 


e p e r f c r rn c . n c c . 


jritten with numerous 


on combustor and no: 
user options to alio w consideration of a number of potential 
combustor and nozzle designs, including provision for jr uel injec- 


tion in liquid f r o rn 1 1 e w c» lie 
pie fuel injector locations ma; 


■(Stream fuel in 


Mu: 


W hi e r "i A D A P T is i mp 1 a n i c n t e d ^ i t ! ■ 1 
raced u r e o r e s u mm a r iced as foil c w a : 


s o 1 u t i o n 


.-•0 evaluating initial profile data for the dependent variables 
for which data have been input . Conipute tlement mass 
fractions from the species mass fractions input- Extract 

from memory the necessary t he r mod y na m i c data for each of 

the species included. In general, Initialize the 

ealeu lation . 


b) 

c) 

-!) 


perform equilibrium calculation at given enthalpy, pressure, 
and ole rn e n t rna s s fractions to obtain species. 


obtain all other p r o f i 1 e information, o . g . clcnsi 1 
h e a t t p h y s i c a 1 coordi n a t e s , v i s c o s i t y , e 1 c . , from 


SpSGxf.G 

1 c o ci e s . 


: he 


change in f h e 


d etc r rn i no 
c o n v ectior 
solut i ori of the go v c rning equatio r 




u s l r. g 


■. n \ e a i « 
ode r . 


return to step (b). 


3. Validation of Present Equilibrium Chemistry Model 


T cj \ a 1 i cJ a t e t h e p r c c e n t equilibrium c h e m _ : t 1 ■ y m o del, 1 1 i e 

r cint chemistry results of ADAPT are compared with TEC results. 


(CoHq) is selected as a test fuel. 1- 
_ The 

selected condition is stoichiometric. The numerical “tsulto of 
ADAPT and CEC ,re summarioed in Table 1 and Table 2 . The flame 
temperature end the mole fraction of the dominant s peel vCCq, 

H o 0 ) are t h e U. e y -ori a b 1 c s •: o c v 1 u o. t c . - s v h o w f 


. -L t- _ 

rpr 


b 1 e f 1 


omami 

OF pool*' QVUYY 



c. 


overall agi .ment for" flame temperatu and species mole 


fractions 

are satis f 

actor y 

r 

. 1 

n P/H cas 

e listed in 

Table 1, 

the 

d e v i a t i o n a 

i 

n 1 1 a me 

t e rn p c r 

atu 

r e _ nd species me:- 1 e 

f r 3 c t ions 

■n r e 

within 0.3 

o. 

0 

and 0.3 

n _ 

0 , ' U _ 

P- - 

ti- ely. I 

n T / P case 1 

isted 1 n 

Tabic 

2 , the dev 

i a t 

Ions of 

-• pt?C l a 

rn 

ole f; act 

i o n: a re w i t 

h i r. 0.2 E 

T „ 

. . 1 1 

terms cf CPU 

time 11, 

ted in 

X . 

i o 

b 1 e 1 an d 

2, ADAPT is 

compu t r 

t i o n - 

ally more 

ef 

f icient 

than C 

EC . 

Especial 

ly for T/P c: 

as a, ADAPT is 

roughly four 

times 

fast: 

r t 

ban CEC. 

These n u rn e 

r ical r c 

suits 

indicate that 

ADAPT c 

art ban 

d 1 e 

1 1 1 e c o m j: 

1 e v equili b r 

i u rn c hi e m 

it try 

problems w 

: t h 

the 1 Ini 

ited p. 

a r t 

lei pc t i rag 

c Kerri 1 c a 1 s pc c i e : a nd 

the 

c onipu t a t x o 

no 1 

e f f i c i e 

ncy . 







4. Conclusions 


N u ru erical results of t h e p r os:nt equ i ] 1 br ium c h o mist r y code 
}a a v e a a a tis F actor y n g r o e n i r '■ t i t S'. CEC r c s u 1 ?: s . These n u m er i c a 1 

results Lndi.ce. te that the present equilibrium solver can handle 
the com p 1 e ■■ ? q k j i 1 i b r c u rn c ! i o n ; 1 . . 1 r y p r o b 1 e r n s uj i t h t ! i e 1 i rn i t e d 

participating chemical '.pec ice end the computational efficiency . 
The fut u r e w o r k c i nc 1 u d e t h c i m pic m e n t o t i o n o f A D A P T equil ibr lu rn 
code in CFD codes such no FPNS a nd MAST. 
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Table 1. 

temperature 


Compa r i s on 
and species 


of CPU tim«_ , predicted flame 
mole fraction for P/H condition 


Present 

Code 


CRAY -XMP Flame 

CPU time (.sec) t em p . ( °K ) 


2213.8 


Mole Fraction 
h 2 o CCU 

0.14298 0.10057 


CEC 

Code 


0.070 


2218.0 


0.14295 0.10037 


Table 2. 

temperature 


Comparison of 
and species mole 


CPU time, predicted flame 
fraction tor f/P condition 


Present 
Cod e 


CRAY -XMP Flame 

CPU time (secj temp.(°K) 


Mole Fraction 

HoO CO, 


0.016 


2218.0 


0 . 14299 0 . 10058 


CEC 

Code 


0 . 0 G 6 


2218 . 0 


0 . 14296 


0 . 10038 



